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<^ ■ Abstract 

Flavor violating processes in the quark and lepton sectors are investigated within 

. a realistic supersymmetric SO(10) x A4 grand unification model. By employing exotic 

heavy fermion fields, this model successfully describes various features of the fermion 

masses and mixings including large neutrino mixings accompanied by small quark 

mixings. In this model the flavor violation is induced at GUT scale, at which A4 flavor 

symmetry is broken, as a consequence of the large mixings of the light fermion fields 

q_) ' with these exotic heavy fields. The stringent experimental constraint from \i — > 

decay rate necessitates a high degree of degeneracy of the supersymmetry breaking 

soft scalar masses of the exotic heavy fields and supersymmetric scalar partners of 

the light fermion fields. The choice of slepton masses of order 1 TeV is found to be 

■ consistent with the constraints from branching ratio of /x — > and with all other flavor 

t^- . changing neutral current processes being sufficiently suppressed. 

O ■ 

1 Introduction 

Flavor changing neutral current (FCNC) processes impose severe constraints on the soft su- 
persymmetric breaking (SSB) sector of the minimal supersymmetric standard model (MSSM). 
The simplest way to satisfy the FCNC constraints is to adopt universality in the scalar masses 
at a high energy scale where the effects of supersymmetry (SUSY) breaking in the hidden 
sector is communicated to the scalar masses of MSSM via gravitational interactions. For 
example, in the the minimal supergravity model (mSUGRA) pQ the MSSM is a valid sym- 
metry between the weak scale and grand unification scale (Mgtjt) at which the universality 
conditions are assumed to hold. In this case, the leptonic flavor violation (LFV) is not in- 
duced. However, in a different class of models studied in Refs [21 [3j HJ O El [7] the universality 
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of the scalar masses will be broken by radiative corrections. Consequently, FCNC will be 
induced in these models as discussed below. 

If the universality conditions hold at the grand unification scale Mgut ; the LFV is induced 
below GUT scale by radiative corrections in the MSSM with right-handed neutrino [2j |3j S] 
or SUSY-S'L r (5) [5] models. Unfortunately, it is difficult to predict LFV decay rates in these 
models because the Dirac neutrino Yukawa couplings are arbitrary within MSSM. However, 
in SO (10) GUT model, we can predict the LFV decay rates below the GUT scale because 
the Dirac neutrino couplings are related to the up-type quark Yukawa couplings and are thus 
fixed. 

The FCNC could also be induced above the GUT scale by radiative corrections. It was 
shown that as a consequence of the large top Yukawa coupling at the unification scale, SUSY 
GUTs with universality conditions valid at the scale M*, where Mqut < M* < Mpi anc k, 
predict lepton flavor violating processes with observable rates [61 [7J. The experimental 
search for these processes provides a significant test for supersymmetric grand unification 
theory (SUSY GUT). Both contributions of FCNC that are induced above and below Mqut 
will be studied in our model. 

In this paper, the flavor violation processes for charged lepton and quark sectors are 
investigated in the framework of a realistic SUSY GUT model based on the gauge group 
50(10) and a discrete non-abelian A4 flavor symmetry [Sj. This model is realistic because it 
successfully describes the fermion masses, CKM mixings and neutrino mixing angles. This 
work differs from other studies in several aspects. First, it is different from those based 
on MSSM with right-handed neutrino masses or SUSY SU(5) in the sense that the Dirac 
neutrino Yukawa couplings are determined from the fermion masses and mixing fit of the 
5*0(10) x A4 model. Thus, this model predicts the lepton flavor violation arising from the 
renormalization group (RG) running from Mqut to the right-handed neutrino mass scales. 
Second, it is different from those based on SUSY 50(10) studied in [9] in the sense that the 
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FCNC processes are closely tied to fermion masses and mixings. Finally, in the SO (10) x 
model flavor violation is induced at the GUT scale at which symmetry is broken due to 
large (order one) mixing of the third generation of MSSM fields (^3) with the exotic heavy 
fields (xi, i runs from 1 to 3). This large mixing arises when the A 4 flavor symmetry is broken 
at the GUT scale. This is different from the case where the flavor violation is induced due 
to large top Yukawa coupling at the GUT scale [5J [7] . The reason for introducing the exotic 
heavy fermion fields in our model is to obtain the correct fermion mass relations at the GUT 
scale as we shall see in section 1. The mass scales of these exotic fields range from 10 14 GeV 
to 10 18 GeV depending on the values of the Yukawa couplings and the scale of A 4 flavor 
symmetry breaking. 

In this paper we study flavor violation of the hadronic and leptonic processes by calcu- 
lating the flavor violating scalar fermion mass insertion parameters (Sab)^ = ^4§^ , for 
(A, B) = (L,R), with m being the average mass of the relevant scalar partner of standard 
model fermions (sfermions). All the flavor violation sources are included in our calculations. 
The sfermion mass insertions, 5ll,rr,lr, arise from the large mixing between the ^3 and Xi 
and the mass insertions, {b % li) RHN , arise from RG running from Mqut to the right-handed 
neutrino mass scales. These scalar mass insertion parameters are analyzed in the framework 
of our model then they are compared with their experimental upper bounds. We found 
that the most stringent constraint on flavor violation comes from the \x — > e-f process. This 
constraint requires a high degree of degeneracy of the soft masses of MSSM fields and the 
exotic fields. Therefore, in this model we assume that these soft masses are universal at 
the scale M* with M* > Mqut-, then we run them down to the GUT scale. The branching 
ratio Br(fi — » e'j) close to experimental bound (i.e. Br(/i — > ej)—1.2 x 10 -11 ) is obtained 
when the slepton masses of order 1 TeV , while the Yukawa couplings remain perturbative 
at the scale M*. We also found in the framework of our model that once the constraint 
from Br(fi — > cy) is satisfied, all the FCNC processes will be automatically consistent with 
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experiments. 

This paper is organized as follows. In section 1, we show how the fermion mass matrices 
are constructed in 50(10) x v4 4 model. In section 2, we discuss the sources of flavor violation 
by finding the sfermion mass insertion parameters 5 l [ L RR at the GUT scale at which A 4 
symmetry is assumed to be broken as well as below the GUT scale. The results of the 
50(10) x A 4 model regarding flavor violation analysis are presented in section 4. Section 5 
has our conclusion. The derivation of the light fermion mass matrices and the light neutrino 
mass matrix after disentangling the exotic fermions is shown in appendix A. In appendix 
B, we list the renormalization group equations (RGEs) for various SUSY preserving and 
breaking parameters between Mqut and M* relevant for FCNC analysis. 

2 A Brief Review of Minimal 50(10) x A A SUSY GUT 

In the 50(10) gauge group, all the quarks and leptons of the SM are naturally accommodated 
within a 16- dimensional irreducible representation. However, minimal 50(10) (i.e., with only 
one 10-dimensional Higgs representation) leads to fermion mass relations at the GUT scale, 
such as ^% = and m° = m®, that are inconsistent with experiment. This can be fixed by 
introducing exotic 16 + 16 fermions and by coupling 16j with these exotic fields via 45//, which 
is used for 50(10) symmetry breaking. The non-abelian discrete A 4 symmetry is chosen in 
our model because it is the smallest group that has a 3-dimensional representation, so the 
three generations of SM fields transform as triplet under A*. Besides, FCNC is not induced 
in the SUSY-5O(10) x A 4 as long as A4 symmetry is preserved. However, as we will see 
later, the breaking of A4 symmetry at the GUT scale will reintroduce the FCNC via large 
mixing between the exotic and light fields. Based on the above reasons, a 50(10) x A 4 
model is proposed in [8]. In this model, a minimal set of Higgs representations are used to 
break the SO(10) gauge group to SM gauge group so the unified gauge coupling remains 
perturbative all the way to the Planck scale. Employing this minimal Higgs representation 
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Table 1: The transformation of the matter fields under SO (10) x and Z2 x Z4 x Z2. 

and symmetry, our model successfully accommodates small mixings of the quark sector 
and large mixings of the neutrino sector in the unified framework as shown summarized 
below. 

The fermion mass matrices of the model proposed in [5] were constructed approximately. 
In this section, we construct these matrices by doing the algebra exactly and show that the 
excellent fit for fermion masses and mixings is obtained by slightly modifying the numerical 
values of the input parameters of Ref.jS]. There are two superpotentials of the model. The 
first one (W spin ) describes the couplings of the standard model fields (^(16^), i runs from 
1-3) with the exotic heavy spinor-antispinor fields (xi(lGi), Xj(16j), % runs from 1 to 3), while 
the second one (W vect .) describes the couplings of ipi with the exotic 10-vector fields (fa, fa, 
fa-, fa", % runs from 1 to 3) as given below: 

W spin , = b 1 4' i XilHi + b 2 ip i X 2 l'm + hXlX3^H + aX3X2lO H + M a XaX a i (1) 

W vect . = b 3 ipi fal6 H + M^fafat + h'^fafajlffk + hijkfafalffk 

+A ijk fa i fa;i'l Ik + mfafa' + hfaffa^H. (2) 

The above superpotentials are invariant under A4 and the additional symmetry Z 2 x Z^x 
Z 2 . The transformations of the matter fields (i.e., the ordinary and exotic fermion fields) 
and the Higgs fields under the assigned symmetry are given in Table [1] and [2j 

The general fermion mass matrix structure that results from integrating out the exotic 
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Table 2: The transformation of the Higgs fields under SO(10)xAt and Z2 x Z4 x Z2. 
heavy spinor-antispinor fields in W spin , is: 

M F (spin.) = f^i?M) f Q 



(3) 



where we have made the following transformation: V^i+V^^+V^ — ^3 an d V'i s i+'02S2 + 
■03 s 3 = 5'('02 s ^+'03 c e)- Here e» and Sj are VEV-components of (Iff) and (1#) respectively and 
s e (c e ) is sin %os 9). f = (1 + T| + T 2 (l + sgTf ))"V2 and r F = (1 + Q 2 F T*T?(1 + s^T 2 )/ 2 ) 1 / 2 
are factors that come from doing the algebra exactly (see appendix A). Here 7\ = 
^2 = , T 3 = ^ and Q = 2I 3 r + §5(t) is the unbroken charge that results from breaking 
SO (10) to the SM gauge group by giving a VEV to 45 H , where (45jj) = fiQ- The charge Q 
for different quarks and leptons is given as. 



1 4 2 

Qu = Qd= gO, Q u c = -1 - -d, g dc = 1 + -5, 

3 6 
Qi = Qn = —-5, Qic = 1 + -5, = -1. 



(4) 



The above general structure of fermion mass matrix has the following interesting features: 
(1) The relation m° = automatically follows from Qd + Qd c — Q e + Qe c , (2) The hierarchy 
of the the second and third masses generation is obtained by taking the limit s e — > 0, and 
(3) The approximate Georgi-Jarlskog relation m° = 3m° s leads to two possible values for 
5, either 5 — > or 5 — >■ —1.25, (4) the former possibility is excluded by experiment since 
it leads to (m^./m^)/(m° s /m^) — > 1 at the GUT scale, while the latter possibility leads to 
(m®/mt)/(m° s /m®) — > which is closer to experiments. Let us define 5 = l + a. The masses 
and mixings of the first generation arise from W vector . The full mass matrices arising from 
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where the parameters are defined in terms of the Yukawa couplings of the superpotential 
(W sp in. + W vect ) and the VEVs of the Higgs fields as shown in appendix A. These matrices 
are multiplied by left-handed fermions on the right and right-handed fermions on the left. A 
doubly lopsided structure for the charged lepton and down quark mass matrices of Eq.(jS]) can 
be obtained by going to the limit (3, (, a, 5s, cu, s < 1 and 6i, 82 are of order one. This doubly 
lopsided form leads simultaneously to large neutrino mixing angles and to small quark mixing 
angles. Based only on the above fermion mass matrices in Eq.Q, an excellent fit is found for 
fermion masses (except for the neutrino masses), quark mixing angles and neutrino mixing 
angles (except the atmospheric angle) by giving the input parameters, appearing in Eq.flS]), 



the following numerical values: 8\ 
s = 0.2897, C = 0.0126, c 12 = -O.OOlle 



■1.28, 6- 



1.124i 



, "2 



and j3 



1.01, 5 3 = 0.015 x e 



1.95/ 



a 



-0.0668, 



-0.11218. The above numerical values 



lead to sin = 0.92 which is not close to the experimental central value of atmospheric 
angle is sin #23™ = 0.707 [10] • This contribution to the atmospheric angle is only from the 
charged lepton sector. Therefore, the neutrino sector should be included by considering the 
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following superpotential: 

W N = b^iZiU H + h ljk Z { Z^ k + rmZfZ?, (6) 

where two fermion singlets Zi and Zf that couple with the Higgs singlet l"^ k have been 
introduced. 

The full neutrino mass matrix is constructed in Appendix B. The Higgs singlet l"^ k has 
the VEV- components (oti, a2, 03). The light neutrino mass matrix is obtained by employing 
the see-saw mechanism. The numerical values (a± = 0.075, ct2 = 0.07, a% = 0.9, and 
A = 0.0465 eV), where A is defined in appendix B, lead to not only the correct contribution 
to the atmospheric angles (sin#23 m — 0.811) but also to the correct light neutrino mass 
differences. The predictions of the fermion masses and mixings are slightly altered by doing 
the algebra exactly compared to the analysis of Ref . [8] . These predictions and their updated 
experimental values obtained from [10] are shown in Table |3j The right handed-neutrino 
masses arise from integrating out the exotic fermion singlets Zi and Zf in Eq.(]6]). The right 
handed-neutrino mass matrix is 



M R = A 



I a$ aitt2l^T + a 2 +a 2 +a 2> a3(a 2 +Q 2 +a 2 ) \ 

n n (-1 1 2 \ ' a| " -azi-al+a'i+al) 



(7) 



where A = 8.45 x 10 15 GeV and the right-handed neutrino masses are given by Mm pa 
M m pa 1.4 x 10 12 GeV and M m = 8.5 x 10 15 GeV. 

Another interesting feature of this model is that it contains a minimal set of Higgs fields 
needed to break 50(10) to the SM gauge group. Consequently, the unified gauge coupling 
remains perturbative all the way up to the Planck scale. This can be understood from the 
running of the unified gauge coupling with energy scale /i > Mgut as 
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Table 3: The fermion masses and mixings and their experimental values. The fermion 
masses, except the neutrino masses, are in GeV. 

where a = g 2 /(47r) and be = S(R) — 3C(G). Here C(G) is the quadratic Casimir invariant 
and S(R) is the Dynkin index summed over all chiral multiplets of the model. The unified 
gauge coupling stays perturbative at the Planck scale (i.e g(M P ) < \/2) as long as be < 
26. Employing large Higgs representations might lead to be > 26. For example, using 
126#+126# gives be = 46. On the other hand, the SO (10) x At-model gives be = 19 which 
is consistent with the unified gauge coupling being perturbative till the Planck scale. 

We will use the same fit for fermion masses and mixings to calculate the mass insertion 
parameters 5 % [ L RR , and 5 % [ R RL in the quark and lepton sectors and consequently investigate 
the FCNC in this model. The charged lepton and down quark mass matrices in Eq.(|5]) are 
diagonalized at the GUT scale by bi-unitary transformation: 

<T = V R ' l M DjL V^' 1 , (9) 

where V R ^ 1 are known numerically. Now, we discuss the sources of FCNC in this model. 
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3 Sources of Flavor Violation in 50(10) x A4 Model 

We assume in our flavor violation analysis that A 4 flavor symmetry is preserved above GUT 
scale and it is only broken at GUT scale. In this case flavor violation is induced at GUT 
scale where A A symmetry is broken. In this section we discuss the flavor violation induced at 
the GUT scale by studying the sfermion mass insertion parameter 5 % [ LRR and the chirality 
flipping mass insertion (A-terms) parameter 5 % l RRL . We will see that these flavor violation 
sources arise from large mixing of the light fields with the heavy fields. This large mixing 
is due to the breaking of A 4 symmetry. In addition, we discuss the induced flavor violation 
arising below GUT scale through the RG running from Mqut to the right-handed neutrino 
mass scales. 

3.1 The Scalar Mass Insertion Parameters 

Let us assume the soft supersymmetry breaking terms originate at the messenger scale M*, 
where M GUT < M* < M Planck . The quadratic soft mass terms of the matter superfields that 
appear in the superpotential W spin , are 

- C = m^l^i + rh 2 Xt x\xi + ^XiXi- (10) 

The MSSM scalar fermions that reside in ^ transform as triplets under the non-abelian 
A 4 symmetry. Since the A 4 symmetry is intact, they have common mass (m^) at the scale 
M*. On the other hand, the exotic fields each of which transforms as singlet under A 4 
symmetry have different masses (rh^,., m|., % runs 1-3) at the scale M*. 

The MSSM scalars remain degenerate above the GUT scale where the A4 symmetry is 
broken. In order to find the scalar masses in the fermion mass eigenstates, two transfor- 
mations are required. The first transformation is needed to block-diagonalize the fermion 
mass matrix into a light and a heavy blocks as shown in Appendix A. The upper left corner 
represents the 3x3 light fermions mass matrix. The second transformation is the com- 
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plete diagonalization of the light fermion mass matrix. Applying the first transformation 
to the quadratic soft mass terms of Eq.l fTUl) by going to the new orthogonal basis (L 2 , L 3 , 
Hi, H 2 , H 3 ) as defined in appendix A, the quadratic soft mass matrix of the light states is 
transformed as follows: 

mp rhp + 5ml, ( n ) 

where, 

/ \ 

Sml = e , (12) 

e = £TZse(ml -ml), 5= ((£) 2 - l)mj + (^) 2 (rr^T? + m^Tf + m^^Tf ), ^d we 
have safely ignored the terms that contain <C 1. It is obvious that the first two generations 
of the light scalars are almost degenerate because the mixing of the second light generation 
(L 2 ) with the heavy states is proportional to s$ <C 1. On the other hand, since the mixing 
of the third light generation (L 3 ) with the heavy states is of order one, its mass splits from 
those of the first two generations. 

The top Yukawa coupling is given in terms of T 1; T 2 , and T 3 as: 

The numerical values of T\ = 0.0305, T 2 = 2, T3 = 100 and a ~ 1.2 are consistent with the 
top Yukawa coupling at the GUT scale to be of order Xf UT ~ 0.5 and r u ^ to be of order 
one. Plugging these numerical values and sq = 0.0465 into the expressions for e and 5 gives 
us: 

(S d ,S d c,S e ,S e c) = (0.81, 0.87, 0.88, 0.82)(m* -mj), 

(e d ,e d c,e e ,e e c) = (0.061, 0.05, 0.048, 0.06) (mj - mj). (14) 

Here we have dropped m 2 Xl terms because their coefficients are negligible. Also, the RGE 
expressions of fh? and m? are the same (see Eq. (J62"]) ). so we have assumed that rr? = 



™L = mi- 



ll 



The next step is to apply the second transformation by evaluating ' 6rr&V L ' and 
similarly for L —¥ R. The unitary matrices Vff are numerically known from the fitting for 
fermion masses and mixings. So, the mass insertion parameters for charged leptons and 
down quarks are given respectively by 

( 5 LL,RR)ij = (VffimljVfyij/mli. (15) 

The above mass insertion analysis without including the superpotential W vect . is good enough 
because we assumed in our analysis that the mixing of the 10 vector multiplets with the 
ordinary spinor fields is small. 

3.2 The Chirality Flipping Mass Insertion (i-terms) 

The FV processes are also induced from the off-diagonal entries of the chirality flipping mass 
matrix Mm- The chirality flipping soft terms are divided into two parts C sp i n and C vect : 

- C spin = bibt^iXilm + hh^ifai'm + hhXiXs^H 

+aax3X2lO H + GiMiXiXi, (16) 
-C vect = b 3 b 3 i!i(j)il6 H + BwM^fa&i + h\jkKjk& j^-Hk + hjkhjk^i^-Hh 

+A ijk A ij kP i <P' j lm + 9™P'i<?"i + hhP'WitiH. (17) 

The fourth term of Eq. (fl6l) induces the off-diagonal elements of the chirality flipping 
mass matrix, if it is written in terms of the new orthogonal basis defined in Eqs.(l26l). This 
transformation can be represented by 

M 2 RL (spin.) aM F (spin.), (18) 

where Mp(spin.) is defined in Eq.([3]). The entire chirality flipping mass matrix in the new 
orthogonal basis is obtained by including —£ vect . The bi-unitary transformations that block- 
diagonalize the full fermion mass matrix is applied on the entire chirality flipping mass matrix 
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(see Appendix A). Accordingly, the 3x3 quadratic mass matrix (M| fl ) associated with the 
light states is transformed as follows: 

M RL -»■ aM F (spin.) + b 3 M F (vector), (19) 

where M F (vect.) = —mM~ l M' (see Eq. (l3~Tj) ) and we have assumed for simplicity that the 
soft parameters appearing in Eq.l fPTt) are all of the same order. Then, the M FR matrix is 
written in the fermion mass eigenstate basis as: 

M RL ->• V^(aM F (spin.) + b 3 M F (vect.))V L . (20) 

It is straightforward to show that the chirality mass insertion parameters are given by: 

{SRLh = %Kt 9 '^ + {zV ] R M F {spinor)V L ) l31 (21) 
m f 

where M d F ia9 - = V ] R M F V L and z = %f . The induced FV arises only from the second term 
of Eq.(T2lD. 

3.3 Mass Insertion Parameters Induced Below Mqut 

The Dirac neutrino Yukawa couplings (ijv)ij induce flavor violating off-diagonal elements in 
the left-handed slepton mass matrix through the RG running from Mqut to the right-handed 
neutrino mass scales. The RGEs for MSSM with right-handed neutrinos are given in Ref. [3]. 
The right-handed neutrinos Mr ( are determined in the 5*0(10) x A4 model. In this case, the 
induced mass insertion parameters for left-handed sleptons are given by [7], 

(Xl \RHN 3m l + ° 2 V^/yr v ]rs M GVT (99 , 

- ^ Y ^ Y Nh* ln ^> ( 22 ) 

where the matrix Yn is written in the mass eigenstates of charged leptons and right-handed 
neutrinos. The total LL contribution for the charged leptons is given by 

(OS * = (JDS™ + (23) 
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Figure 1: The above graphs show the plot of Log of Br(/i — > c-f) divided by experimental 
bound (1.2 x 10~ n ) versus for two cases I and II with M 1 / 2 =787 GeV, 437 GeV and 175 
GeV. 

4 Results 



In this section, we investigate the flavor violating processes by calculating the mass insertion 
parameters 5ll, 8rr, and Slr,rl, then we compare them with their experimental bounds. 
These bounds in the quark and lepton sectors were obtained by comparing the hadronic 
and leptonic flavor changing processes to their experimental values/limits [HI [12]. Eq. (112j) . 
Eq.flHJ) and Eq. (fT5|) are used to calculate 5ll,rr and Eq. (j2T|) is used to calculate Sl^rl for 
both charged leptons and down quarks. The result of mass insertion calculations and their 
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(0.173, 0.016) 


1.7 x 10~ 2 


Wll 


0.054 a 


0.16 


($23) RR 


0.29 tr 


0.22 


{^23)LR,RL 


(0.875, 0.064) k 


(0.006, 0.0045) 



Table 4: The mass insertion parameters predicted by SO (10) x A4 model and their experi- 
mental upper bounds obtained from [T2] . 

~ 9 ~ 9 

experimental bounds are presented in Table HI In this table, we have defined a = x ~ - 2 ^ 

m ii>i 

and k = zm^r 

The stringent bounds on leptonic 612, #13, and #23 in Table S] come only from the decay 
rates Zj — > Ijj. The experimental bounds on the mass insertion parameters listed in column 
3 were obtained by making a scan of mo and M\ji over the ranges m-o < 380 GeV and 
M1/2 < 160 GeV, where mo and M1/2 are the scalar universal mass and the gaugino mass 
respectively [T2] . 

Glancing at Table HJ we note that the stringent constraint on leptonic flavor violation 
arises from 8 12 which corresponds to the decay rate of fi — > ey. On the other hand, there is 
a weaker constraint that arises from Sf 2 on the quark sector. One can do an arrangement 
such that a — 63 = 200 GeV and rhf = 800 GeV (equivalent to k = 2.6 x 10~ 4 ) so that all the 
chirality flipping mass insertions will be within their experimental bounds. This arrangement 
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is possible if the trilinear soft terms vanish at the scale M*. 

Since the stringent constraint comes from the // — > process, let us discuss the branching 
ratio of this process in more details. In general, the branching ratio of Zj — > ^7 is given by 

BRjJj — > Iff) _ 487r 3 a jj 2 A ij\2\ ( 0A \ 
BR{k T^~) " ~GT {1 J + ' Rl } ' [ } 

We have used the general expressions for the amplitudes A l [ R given by Ref. [H] where the 
contributions from both chargino and the neutralino loops are included. These expressions 
are written in terms of mass insertion parameters. 

The correct suppression of the decay rate r(/i — > ej) requires a high degree of degeneracy 
of the soft mass terms of MSSM fields and the exotic fields. For example, a ~ 0.01, as can 
be seen from Table HJ In order to obtain high degree of degeneracy, let us assume that the 
SSB terms which are generated at the messenger scale M* satisfy the universality boundary 
conditions at the scale M* given by 

m \ = m io H = m i H = m l' H = m o, 

2 = 0, (25) 

where M\ is the gaugino mass of S'O(IO) gauge group. Solving the RGE listed in Appendix 
C with the boundary conditions given by Eq. fT2"5]) determines the value of o . In Table [5] 
we give the branching ratio of the process /i — > predicted by the S'O(IO) x A4 model 
for different choices of the input parameters a, bi, 62, and M1/2 at the GUT scale. 
The experimental searches have put the upper limit on the branching ratio of /i — > ey as 
Br(/i — > e^f) < 1.2 x 10 -11 [13J. Note that and Mi/ 2 originate respectively from m and 
Mr, through RGEs. In this Table we consider In -r^- — = 1 and In — = 4.6 that correspond 
respectively to M* w 3M GUT and M* w M P i anck . 

Let us analyze the four cases in the Table El In the cases (I, II and III), the chosen 
values of the parameters a are consistent with the top Yukawa coupling of order 0.5 at the 
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= Ki 




= M , 


a 


= b 1 = 





i 


II 


III 


IV 


a 


1.14 


1.07 


1.14 


0.62 


bi 


1.9 


1.5 


1.24 


1.24 


b 2 


1.9 


1.5 


1.24 


1.24 


Thin 


542 


886 


2932 


675 


M l/2 


350 


787 


1924 


350 


bro -> &y) 


1.4 x 1(T 13 


1.16 x KT 11 


1.2 x 10~ n 


2.2 x 10~ 12 



Table 5: Branching ratio of jj — > cy for different choices of input parameters at the GUT scale. 
Cases I and II correspond to In * f * = 1 and cases III and IV correspond to In ,,f f * = 4.6. 
rhfa and My 2 are given in GeV 

GUT scale and with the fitting for fermion masses and mixing. On the other hand, the 
choice of a = 0.68 in Case IV is not consistent with the fit. Although the medium slepton 
masses of order 550 GeV are obtained in Case I, the choice b\ = b 2 = 1.9 corresponds to 
non-perturbative Yukawa couplings at the scale M* (i.e. b\ = b 2 = 4 at M*). In this case, 
the solutions of the 1-loop RGEs are not trusted since the Yukawa couplings b\ and b 2 go 
non-perturbative above the GUT scale. Also, it is important to point out that the flavor 
violation constraint on fi — > ej in Case III requires heavy slepton masses (> 3 TeV) while it 
requires slepton masses of order ~ 900 GeV in Case II. In other words, Case II is preferred 
in our model in the sense that the decay rate of // — > is close to the experimental limit 
with a reasonable supersymmetric mass spectrum, so it might be tested in the ongoing MEG 
experiment |15j. Besides, the Yukawa couplings remain perturbative at the messenger scale 
M*. Figure [T] shows the allowed values of that correspond to the graphs below the x-axis 
for the cases I and II. 

5 Conclusion 

In this paper, we investigated flavor violating processes that arise below and above the GUT 
scale in the 5O(10) x model. Above the GUT scale, we study how flavor violation gets 
linked with the fitting of fermion masses and mixing through the factors T 1; T 2 , and T 3 . 
The requirement of top Yukawa coupling being ~ 0.5 necessitates some of these factors to 
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be large. Consequently, this corresponds to an order one mixing of the light fields with the 
exotic heavy fields. In this case, flavor violation is reintroduced at the GUT scale where 
A4 symmetry is broken. The stringent constraint on fi — > ej decay rate requires a high 
degree of degeneracy of the soft quadratic masses of the exotic heavy fields and the light 
fields. Therefore, all the quadratic soft masses are assumed to be universal at the scale 
M* ~ ?>Mgut- Flavor violation is also induced below the GUT scale in the presence of 
right handed neutrinos through the RG running from Mqut to the right handed neutrino 
mass scales. This FV source is predicted by SO(10) x A4 model because the Dirac neutrino 
Yukawa couplings are determined from the fermion masses and mixings fitting. Combining 
all sources of FV, we found that Case I and Case II presented in Table [5] are consistent with 
fermion masses and mixing fitting and with \i — > cy decay rate, which is however predicted 
to be close to the current experimental bound. Thus the ongoing MEG experiment can 
confirm or rule out our model. Case I that corresponds to slepton mass of order 1 TeV is 
also consistent with the Yukawa couplings (i.e., b\ and b 2 ) being perturbative at the scale 
M*. On the other hand, these Yukawa couplings do not remain perturbative at M* in Case 
II that corresponds to slepton masses of order 550 GeV. 
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6 Appendices 

A Derivation of the Light Fermion Mass Matrix 



In order to block-diagonalize the mass matrix of W sp i nor} we define the new orthogonal basis 
as Y — UX, where Y(X) is the column matrix that contains the new(old) eigenstates and 
U is the 5x5 orthogonal matrix (i.e U T U = UU T = I). These matrices are given by: 



/ L 2 \ 




La 




#1 













\ 



-Ni 








N lSe T 2 





fN lCe s e T^ 


/ 




fN lCB T 2 




r F 




Nir F 


r F 


Ntrp 








GfQfTz 





Gp 


N 2 sgT 2 


N 2 CgT 2 





^2 





fN 2 c e s g TiT 1 


fN 2 T ± 


fG P 


fN 2 c e T 2 T t 


fQpG F T 3 


G f tf 


N'{G F r F 


N 2 r F 


Gprp 


N 2 r F 



\ ( ife \ 

Xl 
X2 

) \ X3 / 



(26) 



where N ± = 1/^1 + 1^, N 2 = 1/y/T+T?, G F = l/y/1 + T*Q% f = (1 + If + T 2 (l + 
SgT|)) -1 / 2 , and rjr = + Q 2 F T^T^(1 + SgT 2 2 )/ 2 ). The parameters appearing in the above 
matrix are assumed to be real. Define E iy E c i) g i} g[, g", and g"' to be the charge (— 1) 
leptons in the ip h \u Xi, <f>» <t>'i, <P", and respectively; and define e£, E h gf, g'f, g'/ c , 
and g'l' c to be the charge (+1) antileptons in the same representations. By writing the old 
eigenstates appearing in the superpotential (W spin + W vect ) of Eqs.flT]) and © in terms of 
the new ones, and restricting attention to the electron-type leptons, one gets a 21 x 21 mass 
matrix: 



W„ 



e i ex. h a 9i 9i 9i 9i ) 



mo m 
M' M 



/ 


ei 






E a 












9i 






9i 






9'1 




V 


9? 


1 



(27) 



where, 





(0 





m = 










1° 


afVdQe cS B 







\ 
J 



The matrices M', m and M can be written in the compact form as 
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M' 



m 











M = (Mn M12 M13 ) , 



(28) 

(29) 
(30) 



where 



Ki = 
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V o 
/ o 








b 3 v 5 



V 



Nir e 



aN 1 v d G e cs e T 2 


"He 
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N 2 T e c 
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-6 3 iVit; 5 
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afN 2 VdQ e cT 1 T 3 
Nir e c 

apv d c e {NlGlQ e -N? 2 Q e cT^)T 2 T 3 

G e r e r e c 







bsfNiV5CeseTl 
r e c 

N!r e c 



J 
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aN 2 v d G e c 
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and 



M 13 = 
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k 2 riQ e 








m 








o / 



Here v% = (1(16//)), v 5 = (5(16//)), Vd = (5(10//)), s e = sin6> and ce = cos#. The above 
21 x 21 mass matrix may be block-diagonalized as follows [15] : 



T ( mo to \ TT ] _ ( (to - toM _1 M')(1 + y f y)~ 1/2 \ 

R I M' M ) L ~ \ (M Jlf t + M'Aft) J ' 



where 



T j ( I (TO M' t +TOAft)(AO/t + Af'Af't)- 1 \ 

fl ~~ I (MMt +M'M't)- 1 (TOjM' + TOtM) / y' ^ 



and 



y (MAft + M'M't)- 1 ^ j ^ M' M y ' 

Here y = M~ l M' . Terms of order (Mweak/Mcur) 2 have been dropped. Then the 3x3 
light fermion mass matrix of charged leptons in Eq.(jS]) is obtained by applying the relation 
in the left upper block of the matrix in the Eq. flSTj) . where the factor (1 + y^y)^ 1 ^ 2 is close 
to identity for small mixing between the ipi and the 10-plet vectors. Similarly, one can 
obtain the down- type quark mass matrix. The parameters appearing in Eq.(j5]) are defined 
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as follows: 



C = c l3 + 5 2 ^^, (34) 
= c 23 -5 1 ^^, (35) 



5s e 



f(2 + 3a)c e '' 



(36) 



° 12 af 2 v d c e Mf Q (Q e + Q e c)TiT 2 T 3 ' 1 j 



(^i - A 2 )&|A; 2 Ar 1 t; 1 t; 5 73^ 
af 2 mv d c 9 Mf (Q e + Q e o) 7\T 2 T 3 : 
&jj/i"i"5 (e 2 - Nle z ces 9 Tf: 



^1 - ^ 2y |U 3 A, 2 iV ii/ic/573" /oo\ 
d 3 - r7xZ7 I ~~2 T7\ ~ \rrrrm , 



13 afNiVtCeM^iQe + Q^TtTzTs' K } 

s (Ai - A 2 )b 2 3 k 2 v l v 5 Q (72 + Nf 73 c e s e T^) 



-6 3 /ifif 5 ei 
afv d c M* Q (Q e + g e c)TiT 2 T 3 



tfi = (-Ax + A^blk^v^Q 

afmv d c e Ml (Q e + Q e c) T ± T 2 T 3 

The above parameters are written in terms of the Yukawa couplings and the VEVs of the 
Higgs fields appearing in the superpotentials W spin and W vect . in Eqs.flTJ and (j2J). The 
parameters 71, 7 2 and 7 3 appearing in the above Eqs. (j3"gj) . (jiDj) and (jl2"j) are the VEV 
components of the Higgs singlet 1" H . 

B Light Neutrino Mass Matrix 

The neutrino mass matrix can be obtained from the superpotentials given by Eqs.dTJ and 
OH]). For simplicity, the contribution from the superpotential W vect . in Eq.fll - ) is ignored by 
assuming the coupling of the ordinary spinor fields 16j with the vector multiplets is small. 
Define the right- and left-handed neutrinos, denoted respectively by (uf and Ui), residing in 
ipi. Similarly, v'L and v Xi (v c x. and F-J reside in Xi{Xi) where i runs from 1 to 3. Including 
the six singlets denoted by and Z?, one can construct 24 x 24 mass matrix written in the 
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following compact form 



where 
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' 1 Ml M R ) V 
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I 



Here v u = (5(10 #)). The matrix Mr can be written in the compact form 



M R = ( Mflii Mfli 2 M fl i3 M m4 ) , 



where the matrices Mru, Mri 2 , Mri 3 , and M^i 4 are given respectively by 
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(47) 



Here v = (1(16#)). The light neutrino mass matrix is given by the seesaw formula as fellows 

/ \ 

My = M D M^ l Ml = A k 77 , (48) 

V 77 1 / 
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where 



A Aa 2 c 2 f 2 v 2 T 2 T 2 ((a 2 + a 2 ) Q 2 u r 2 uC s 2 + 2N 2 a 2 a 3 c 6 Q u r u cs e ((Q„ + Q v .) r v 



+ 
V = 



m\N 2 v 2 b\r 2 r 2 c 

Q u r u cS 2 T 2 ) + N[ {a 2 + a 2 ) c 2 {{Q u + Q„ e ) r v + Q u r uC s 2 T 2 ) 2 ) T 2 

miNlv 2 h\r 2 v r 2 c 

NlQ u -r 2 u s e (a 2 a 3 Q u r v c Se + N 2 (aj + a 2 ) c e {{Q v + Q u c) r u + Q u r u c S 2 T 2 )) 



f(A + B) 



Here the numerical values of a±, a 2 , a 3 and A are given in section 2, and we have defined 

A = (a\ + af) Q\ r \* s \ + 2N 2 a 2 a 3 c e Q e r e cs e ((Q e + Q e c) r e + Q e r e c.s 2 e T 2 ) , 
B = Nt{al + a 2 2 )c 2 e ({Q e + Q e c)r e + Q e r eC s 2 e T 2 ). 

C RGEs from the scale M* to the GUT scale 

Neglecting all the couplings in the superpotential W vector , since they do not contribute to 
the top Yukawa coupling, we present only the RGEs that are needed to find the parameter 
a at the GUT scale. The one-loop RGE's of the unified gauge (go) coupling, the couplings 
appearing in W sp i nor , and the trilinear soft terms associated with W sp i nor between the scale 
M* and GUT scale are 

= I9gl (50) 

16tt 2 ^ = 6 1 (206 2 + 6 2 -45^), (51) 

16vr 2 ^ = 6 2 (206 2 + 6 2 -45^), (52) 

16^ = a(l 8a 2 -f, 2 ), (53) 
db 

16tt 2 ^ = 2(206 2 6 1 + 6 2 6 2 + 45^M A ), (54) 

16tt 2 ^ = 2(20b 2 2 b 2 + b 2 1 b 1 + A5g 2 G M x ), (55) 



dt 

16tt 2 ^ = 285a 2 + 63g%M x . (56) 
dt 
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The RGE's soft mass terms for the fields appearing in W sp i nor are given below: 



16vr 2 ^ = 2^(174+7^+7^+6?) 

+ 26*(tt4 + 774 + ml + 6 2 ) - Ahg 2 G Ml (57) 

cItyi— ~ 

l ^ 2 —^ = 66?(m^+m2 i+ m^+6?)-45^M A 2 , (58) 

diri— ~ 

16?r2 ^f = Qb lH,+ rh2 x 2 + rh2 i' Ht +~ b l)- 45 9 2 G Ml (59) 

1Q ^^r = 326 2 (mJ i +m| i +m 2 fl! +6 2 ), (60) 
dml 

l ^ 2 ^T = 32bl(ml+ml 2 +m 2 % +b 2 2 ), (61) 

16n2 ~i £ = 10a2 K 2 +^x 3 +^o fl +« 2 )-45^M 2 , (62) 

167t 2^+0h = 16a 2^2 +fh/ 2 +lfft 2 + a 2 )-36g 2 G M 2 x . (63) 



Here m\ , ml and m 2 are the quadratic soft masses for the Higgs superfields appearing 
in W sp in defined in Eq.fTJJ) and the quadratic soft masses m^., m 2 7 12 , and ra^ia are defined 
in Eq.dEU. 
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